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Abstract: We discuss explicit examples of BPS solutions in four-dimensional N = 2 su- 
pergravity with i? 2 -interactions. We demonstrate how to construct solutions by iteration. 
Generically, the presence of higher-curvature interactions leads to non-static spacetime 
line elements. We comment on the existence of horizons for multi-centered solutions. 



X 

1 Introduction 

This note is a follow-up of a recent study of stationary BPS solutions in four-dimen- 
sional N = 2 supergravity theories with higher-curvature interactions and with an arbi- 
trary number of abelian vector multiplets and neutral hypermultiplets. Working within 
the framework of the superconformal multiplet calculus an analysis of solutions with 
N = 2 and residual N = 1 supersymmetry was presented. It was established that 
there exist only two classes of fully supersymmetric vacua, namely Bertotti-Robinson 
and Minkowski spacetime, and that the corresponding solutions are fully specified by the 
charges carried by the field configurations. As a consequence of the uniqueness of the 
horizon geometry the so-called fixed-point behavior is thus established in the presence 
of i? 2 -terms. Furthermore, concise equations that govern stationary BPS solutions were 
derived. In the absence of i? 2 -interactions our conclusions agree with previous results 
(here we refer in particular to the recent work ||). In this note we demonstrate that the 
construction of BPS solutions in N = 2 supergravity with higher-curvature interactions 
is in general very complicated and show, by studying a few explicit examples, how to 
proceed iteratively. We address the fact that multi-centered BPS solutions are non-static 
due to the presence of i? 2 -interactions, and present arguments why these configurations 
may still have regular horizons at the charge centers. 



2 Residual supersymmetry 



For what follows we adopt the definitions and conventions used in 0. The line element 
of a stationary solution is parametrized by 

ds 2 = -e 29 (dt + a m dx m ) 2 + e~ 29 g mn dx m dx n . (1) 

Imposing residual N = 1 supersymmetry subject to the embedding condition 

h €i = eij 7o e J , (2) 

where h is a phase, leads to stringent restrictions. For instance, g mn needs to be a flat 
three-dimensional metric and the hypermultiplet scalars must be constant. To transpar- 
ently present the restrictions of residual supersymmetry on the vector multiplet sector it 
is advantageous to rescale the vector multiplet scalars and to express the Lagrangian in 
terms of the scale and U(l) invariant variables Y 1 (I — 0, ... ,n), which are shown to be 
expressed in terms of harmonic functions if 7 (x) and Hj(x) according to 

( Y '-V..)=i( H '). (3) 

These 2(n+ 1) equations for the real and imaginary parts of the scalars Y 1 constitute the 
so-called generalized stabilization equations. In the absence of i? 2 -terms these equations 
were first conjectured in Q. Here Fj(Y, T) is the partial derivative with respect to Y 1 of a 
homogeneous holomorphic function F(Y, T). This function characterizes the couplings of 
the vector multiplets. The field T appearing in F(Y, T) as an extra holomorphic parameter 
is the (rescaled) scalar of a chiral background superfield which, upon identification with the 
square of the Weyl superfield, parametrizes the i? 2 -interactions. The harmonic functions 
are related to the electric and magnetic charges carried by the field configurations. In the 
case of multiple centers they are superpositions of the harmonic functions associated with 
the electric (<7at) an d magnetic (p T A ) charges carried by the centers located at xa, 

^ = ^+EiAt- ^^i+Ef^m > ( 4 ) 

where the (h 1 , hj) are constants. It is remarkable that apart from @ the only other 
equations that must be solved to fully specify the stationary BPS solutions, are the ones 
for the spacetime line element, 

z [y J F 7 (y,T)-F/(y,f)y J ] +|xe- 2s = 128ie 9 V p [(V p e~ 9 ) (F T - F?) 

-32*e 4 M#(a) p ) 2 (F T -F T ) 
-64e 2!? J R(a) p V p (F T + 

Ft) i 

H'VpHj+ixR^p = -128tV[v [p (e 29 R(a) q] (F r -F r ))' 

-128V 9 [2V^V 9] (F T + F T )], (5) 

where T is constrained to take the form 

r = -G4(V p g- 1 1 ie 29 R(a) p ) 2 . (6) 



Here p, q denote three-dimensional tangent space indices and R{(j) p is the dual of the 
field strength of a, R(cr) p = e pqs V q a s . Furthermore, Fx = d~[F{Y, T) and x denotes a 
negative constant which fixes the scale. All other fields, such as the field strengths of the 
vector multiplets, are given by explicit formulae in terms of g, R(a) p , and the moduli 
fields Y 1 . Solutions to the above equations automatically satisfy all equations of motion. 



3 Examples of BPS solutions 

The equations (|]), fl5|), and (j^j that determine stationary BPS solutions can be solved 
explicitly only in very few cases. For instance, it is often difficult to obtain the solution to 
the generalized stabilization equations @ in closed form, even in the absence of i? 2 -terms. 
Furthermore, equations (||) and @ are coupled and can usually be solved only iteratively. 
For concreteness, let us consider a first simple example, 

F(Y,T) = -\iY I r llJ Y J + cT, (7) 

where i] is a real symmetric matrix and c a complex number. In this case it is simple to 
solve the equations (|3|), 

Y 1 = i (iH 1 - V IJ Hj) , Fi = \ (iHj + Vl jH^ 

where f]ijr] JK = 5j K . The dependence on the i? 2 -background will enter only when solving 
for the line element, as we will discuss shortly. This is rather the exception than the 
rule. Consider, for instance, the coupling function which arises in Calabi-Yau threefold 
compactifications in the large- volume limit, 

D Y A Y B Y c Y A 
F(Y,T)= ABC y0 + d A — T, (9) 

with A,B,C = 1,2, ... ,n. We construct solutions to this model satisfying H° = so 
that Y° is real. Introducing the matrix Dab = DabcH c and assuming its invertibility 
DabD bc = &a , the stabilization equations can be solved to all orders in T, 

y a = i Y °D AB (H B + id B T 




(Y o )2 _ DabcH a H b H c - i(T - Tfd A D AB d B - 2(T + T)d A H' 

4 (h + ±H A D**H B ) 



For more complicated F(Y, T), solving (|3|) can become very involved. It may be possible 
to cast F(Y, T) into a power series expansion (possibly after an electric/magnetic duality 
transformation) in which case the generalized stabilization equations (§) can be solved 
iteratively in powers of T, 

Y\H, T, T) = J2 Y Lm)(H)r n t m . (11) 

n,m 

Clearly, Fj(Y, T) and Fy(Y, T) will have corresponding expansions in T and T once the 
solutions Y^H, T, T) are inserted, and one could in principle, by treating T as a formal 
expansion parameter, proceed to solve (§) and (H). Since such a procedure is not feasible 
in practice, the question arises whether it makes sense to solve the equations (||) and (|6]) 



iteratively and to truncate at some suitable order. To address this question we recall that 
the function F(Y, T) is homogeneous of degree two, 



F(XY, X 2 T) = A 2 F(Y, T) . (12) 

It follows that 

F 7 Fr(y, T) + 2T F T (Y, T) = 2F(Y } T) , (13) 

and in particular we have Fi(\Y, A 2 Y) = XFj(Y, T). This shows that the equations 
(§) are invariant under this rescaling if we let the harmonic functions H 1 and Hj scale 
with weight one. Therefore the coefficient functions Y/ nm JH) in the expansion ([TT]) will 
scale with weight 1 — 2(n + m), such that every power of T is accompanied by a net 
amount of two inverse powers of harmonic functions H~ 2 . The expressions ( |i~0f ) illustrate 
this feature. In a similar way homogeneity is reflected at the level of the Lagrangian. 
Therefore, corrections due to i? 2 -interactions become subleading whenever |T| <C H 2 . 
One can pinpoint such a hierarchy when one encounters supersymmetry enhancement in 
the neighborhood of a charged center. There, |T| has a l/r 2 -fall-off proportional to a 
charge-independent constant, while the harmonic functions fall off as Q/r, where Q is the 
charge carried by the center. This is the reason why the corrections to the entropy of BPS 
black holes 0, for instance, are subleading in the limit of large charges. In fact, owing 
to homogeneity, the entropy will have an expansion of the form S = vr £ n>0 S^Q 2 ~ 2n , 
where the coefficients S^- 71 ' are independent of the charges. According to the arguments 
presented above, homogeneity implies that the expansion of the line element takes the 
schematic form 

e" 29 ~ £a (n) |rr# 2 - 2 ", R{a) ~ ]T P (n) \T\ n H 2 ~ 2n , (14) 

n>0 n>0 

where ar n ) an d are independent of the harmonic functions. This shows that a trun- 
cation at some finite order may be sensible only in situations in which T falls off much 
stronger than the harmonic functions if 2 when moving away from the centers. 

Let us reconsider the holomorphic coupling function (0). In this simple example we can 
use c as an expansion parameter, since by homogeneity every power of c will always be 
accompanied by two inverse powers of harmonic functions. After solving the generalized 
stabilization equations @ the remaining equations @ reduce to 

(H^uH 3 + H I r ] IJ Hj)+ X e- 29 = 256 i (c - c) [e 9 V 2 p e~ 9 - (|e 2 ^(a) p ) 2 ] , 

2H 1 V P Hi + X R(cr) P = -256 z(c-c) V q [V [p (e 29 R(a) q] )} . (15) 

The case where c is real corresponds to adding a total derivative term to the action. Above 
formulae show that the line element stays unaltered in this case, while T is given by 



T = 64 



H 1 - t V IJ Hj) V P H! - (Hj + tVijH J ) W p H r 2 
WruJIF + H lV IJ Hj 



(16) 



There are other less obvious situations where the dependency on c drops out of the 
equations fllsp . This is the case, for instance, when R(a) p = and all the harmonic 
functions are proportional to e~ 9 , 



H^^e' 9 , H I = a I e~ 9 , (17) 



where (a 7 , a/) are constants. 

Generically, however, the right-hand sides of dl5|) do not vanish and one has to rely 

on an iterative approach. For mutually local charges, H 1 V P Hi = 0, static solutions 
with R(a) p = are possible. This is what we want to investigate in the following. The 
remaining equation, 



e- 29 (x - 256 i(c - c) e 3ff Vje" 



is a non-linear differential equation for e 9 . To zeroth order in c we find 



-291(0) 



-X 



-i 



(19) 



Making the ansatz e 9 = J2n>oi e 9 Y n \256i(c — c)/x) n the line element is determined 
iteratively by 



-9]( n ) 



\ [e 29 ] {0) (v 2 [e- 9 ] (n -V - £ V 9 ] {^) [e~ 9 ] V 



^7l« 



(20) 



where the truncated sum ij t k runs over a ^ < i, j, k < n subject to % + j + k = n. The 
presence of the overall factor [e 2s ] ^ ~ H~ 2 on the right-hand side of ( p0|) indeed induces 
the expansion indicated by (|14D. 

Let us return to the more complicated example (|9]). In this case we can use the d A as 
expansion parameters. The exact solution to the generalized stabilization equations for 
the case H° = are given in fllPf). We find by direct calculation that 



F 7 F 7 - r 7 F 7 



D ABC H A H B H C 



Y° 



+ d A H- 



T + T 



and 



Fy — Fy = - — — — — , F~f + Fy = i d A D AB \ H B -\- ids 



yo ' 



yO ' 1 i 3 

To zeroth order in d A the solutions to (H) are given by 
D ABC H A H B H C 



[yo](o) 



where [y°]W is the zeroth order expression for Y°, 

D ABC H A H B H C 



([Y 



0l(0)\2 



4(# + ±H A D AB H B ) ■ 



(21) 



:22) 



(23) 



(24) 



Keeping track of the terms coming from ( p^) the expressions for the line element up to 
second order in d A are readily obtained from (W), 



D ABC H A H B H C 



[Y 



01(0) 



1287 _1 V P 



(V„7) 



d A H A 



[yo](o) 



+ 32 7 - 4 ( P} 



i 2 d A H A 



[yo](o) 



fr 2 P P V p (d A D AB H B )+0(d A d B ) 



H 1 V P Hj + 



256. 



7 _1 (V [p 7)V g] (^ J D AS J ff B ) 



+128V 9 V [P 7 "%] 



[yo](o) 



C>(d A d 



(25) 



where j 2 and p p are abbreviations for the zeroth order results as given in fl23|) for e~ 2g and 
R(a) p , respectively. We see again that the first-order approximation to the line element 
is cast into the form (|T^), the first correction from the i? 2 -terms being suppressed by one 
power of H~ 2 . 



F(Y,r) = -^(Yy + b— -. (26) 



4 Static versus stationary BPS solutions 

As pointed out in JIJ, higher-curvature interactions induce non-static pieces in the line 

element even for BPS configurations with mutually local charges, H 1 Vp Hj = 0. This 
effect can be observed in the previously discussed example, but it does not arise when the 
solution has only one center because the right-hand side of the second equation (||) vanishes 
due to rotational symmetry. Therefore the question arises whether i? 2 -interactions still 
allow for regular multi-centered black hole solutions. Again, it is difficult to address this 
question in full generality. Due to fll4]) we expect e~ 2g generically to diverge cLS J X X A | 
as one approaches one of the charge centers xa- When the charges are not mutually local 
one expects R(cr) p to behave cLS X X A 

| . For mutually local charges, on the other 
hand, the singularity of R(cr) p is, in fact, milder. Inspection of the expressions for the 
curvature components given in ffl] shows that every term involving R(a) p is accompanied 
by a sufficient amount of factors of e 9 such that the effect of the non-static pieces of the 
line element disappears near the centers. This indicates that multi-centered BPS solutions 
still possess an AdS2 x S 2 geometry there. The following example illustrates this. 

Let us consider a simple model describing pure supergravity with the particular R 2 - 
interactions given by 

We assume b to be a real constant. By the same arguments as given in the previous section, 
we use b as an expansion parameter. We solve the generalized stabilization equations for 
the purely electric situation H° = 0, H = H. To zeroth order in b we find 

e~ 29 = -x- 1 H 2 + 0(b) , R(a) p = 0(b), T = -Q4H~ 2 (V P H) 2 + 0(b) . (27) 

To calculate the next-order corrections to the line element one needs + Fy to first 
order in b, 

F r + F r = -2 10 b H- 4 (V P H) 2 + 0{b 2 ) . (28) 

We consider a harmonic function H with multiple centers located at xa- For simplicity 
let us assume that center A — 1 is at X\ — and calculate R(cr) p around this center. 
Therefore we expand the harmonic function appearing in above expression in powers of 

V P H = - qi §- 3 + ElA^- 3 +0(\x\), 
V p V q H = Ql {& - + £ 1A ( 3^%^ - r^] + C(\x I) . (29) 

\ \ X \ \ X \ ) Ay^l \ \ X M \ X M J 



Near \x\ ~ one finds that, to first order in b, R(a) p is given by 

R(a) p = ^ £ <U fef + <^ 2 ) , (30) 

AVl H \ X A\ 

where x p denote the components of the unit vector. Thus, R(a) p exhibits only a 
singularity so that the Riemannn curvature near the center A = 1 is not affected. Hence, 
in this simple example we recover the usual AdS2 x S 2 geometry typical for extremal black 
holes. 
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